Formation of helical states in wormlike polymer chains 
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We propose a potential for wormlike polymer chains which can be used to model the low-temperature 
conformational structures. We successfully reproduced helix ground states up to 6.5 helical loops, 
using the multicanonical Monte Carlo simulation method. We demonstrate that the coil-helix transi- 
tion involves four distinct phases: coil(gaslike) , collapsed globular(liquidlike), and two helical phases 
I and II (both solidlike). The helix I phase is characterized by a helical structure with dangling loose 
ends, and the helix II phase corresponds to a near perfect helix ordering in the entire crystallized 
chain. 
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The structure of a-helices is one of the most common 
conformations occurring in biological molecules such as 
proteins. It is this unique structure that acts as the 
building block for much more complicated tertiary struc- 
tures [jp. In 1959, Zimm and Bragg proposed a mean- 
field model for the coil-helix transition that was based 
on the thermodynamics of a known sequence of pairing 
residues A tremendous number of theoretical works 
have been published since then and the basic approaches 
are derivatives or alterations of the original Zimm-Bragg 
treatment |||| . Most recent numerical attempts to study 
the problem are based on more realistic potentials at the 
atomic level for helix- forming polypeptides, concentrat- 
ing more on finding the ground state structure than 
on the thermodynamics at finite temperatures (see, how- 
ever, Ref. [§). 

From the statistical-physics perspective, it is desirable 
to start from a simple theoretical model of general char- 
acteristics, in the same sprite as recent models for double 
stranded semiflexible polymers [Q, for low-temperature 
lattice homopolymers || and for polymer chains con- 
taining dipolar segments |j) . By ignoring the specifics of 
the amino acids, most proteins can be viewed as worm- 
like polymer chains. One can then ask simple questions 
such as, what are the minimum conditions required for 
a wormlike polymer chain to reproduce a helix ground 
state? What are the characteristics of the so-called coil- 
helix transition in comparison with the predictions of a 
Zimm-Bragg-type theory, and how different is it from the 
coil-collapsing transition in other polymeric molecules oc- 
curring near the theta temperature? From a polymer 
physics point of view, a related and inspiring question by 
itself is: what are the low temperature structures of a 
wormlike polymer chain in general? 

There arc two basic characteristics common to all 
wormlike chains: the existence of a finite persistence 
length, which prevents the polymer from folding as a 
sharp hairpin loop at scales smaller than the persistence 
length, and an excluded volume interaction between the 
monomers, which prevents polymer chains from collaps- 



ing into a point at low temperatures. For helix- forming 
wormlike polymers, a third basic characteristic must be 
included: a directionally biased attractive potential. 

In the current approach we model wormlike polymers 
by fixing the bond angle 8 between two adjacent poly- 
mer bonds of length a, and allowing the azimuthal angle 
<j) rotate freely, as prescribed in the original freely ro- 
tating model | flo|| . Movement of monomers is realized 
by implementing a pivot algorithm [pd| , where a bond 
is chosen at random as a reference axis about which the 
entire polymer segment attached to this bond is rotated. 
The Metropolis method of sampling a canonical phase 
space is known to be ineffective for examining the possi- 
ble conformational states at low temperatures, as a par- 
ticular configuration could easily become trapped in a 
local energy minimum. To avoid this problem we have 
incorporated a multicanonical technique into the pivot 
algorithm, which involves reweighting the temperature 
as a function of energy, in order to produce a flat simula- 
tion histogram across the entire energy landscape ||[l2| • 
This approach allows the polymer to potentially visit all 
possible energy states and to tunnel through local energy 
minima. For a more detailed explanation of the multi- 
canonical technique the reader is referred to Refs |],|l2j ■ 

The interaction energy between the monomers labeled 
i and j having a center-to-center distance is proposed 
to have the simple form, 
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for a < Tij 

(n ~ fj)] m - e[uj ■ (fi - rj)] m for d < r z] < a 

for < ra < d 



where Ui is the cross product defined as 

U\ = (n - fi-i) x (r i+ i - ri). 



(1) 



(2) 



Here, m represents the degree of directional bias in the 
potential, d is the excluded-volume diameter, and a is 
the attractive-force range. The multicanonical technique 
is also an effective method to search for the ground states 
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of a given potential. In an effort to answer the question 
of what type of potential is needed for the creation of 
a perfect helical ground state, the consequences of using 
various values of m were examined. We found that for 
an isotropic interaction (to = 0), the stable ground state 
structure is a disordered globular state with no particu- 
lar relative directional ordering between the chain bonds. 
Thus, a directionalized interaction (to ^ 0) is obviously 
an essential condition for a helix-forming wormlikc chain. 
We also found that weaker powers (to = 2,4), would only 
produce a perfect helix ground state in short polymers; 
in longer polymers multiple bonds formed by spooling 
around a helix core i might give rise to an even lower 
energy. In a helix state, bonding effects should be max- 
imized when all monomers siting directly on top of each 
other, demanding a strongly directionalized interaction 
potential as in known helix-forming polypeptides. As an 
illustration of some of the specific physical properties of 
a helix-forming wormlike chain, we report here the case 
of to = 6 only. Numerical data collected on using other 
power to will be reported elsewhere. 

It is worth noting that our potential is not simply 
a function of the distance between the two involved 
monomers. The bonding preference also involves the di- 
rections of the bonds connected to these two monomers, 
which is an interesting effect in the transition to the he- 
lix states, similar to a protypical hydrophobic interac- 
tion in proteins. There are a number of other variables 
in the potential model. The diameter d and constant a 
have been adjusted to have the values 3/2a and y / ^45/8a 1 
while the nearest neighbor interaction is ignored. The 
fixed bond angle 9 is chosen to have a value of approx- 
imately 7r/3, corrected for the helical pitch to maintain 
six monomers per helix loop. These selections of values 
were based on the observation that under the current 
choice a stray end of the monomers cannot loop through 
an open core of a partially formed helical segment. We 
have also investigated the case of = 7r/5, and found no 
major qualitative differenced in the physical properties 
to be examined below. This choice of potential, with the 
use multicanonical technique, was able to produce a per- 
fect helical ground state relatively quickly (see Fig Id) 
for various polymer lengths N = 13, 19, 26 and 39, show- 
ing that the helix ground state can be constructed with 
the above mentioned three conditions of a helix-forming 
wormlike chain. 

Therefore, a wormlike chain can have numerous types 
of ground states, which are dependent on the type of 
directional potential used. The construction a poten- 
tial which produces a singular ground state of a desired 
conformation is possible, although not a trivial exercise. 
The use of multicanonical technique makes the search for 
these specific ground states simpler. Knowing that our 
model has an attainable perfect helix ground state, we 
have calculated the multicanonical reweighting function 
for the temperature range discussed below, which enable 
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FIG. 1. Snapshots of the configurations of a 26mer at var- 
ious temperatures: (a) fcsT/e = 4 (coil), (b) fesT/e = 1.3 
(globular) (c) k B T/e = 0.8 (helix I) and (d) k B T/e = 0.1 (he- 
lix II). The size of the beads represents the actual hard-core 
interaction between n on-a djacent monomers, and the attrac- 
tion force range is y^5/2a, where a is the bond length. A 
highly directionalized potential with m = 6 (see Eq. (1)) is 
used. 



us to conduct measurements of physical quantities in the 
subsequent production runs. As a most effective measure 
for exploring the phase transition signatures, the spe- 
cific heat capacity of the polymer chain was calculated 
by measuring the energy fluctuations 
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The mean square radius of gyration, 
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< R i >=< jf i> - fcm ) 2 > 

i=l 



(3) 



(4) 



where r cm is the center of mass vector, and the mean 
square end to end distance, 



< ^end >=< (n - ^n) 2 > 



(5) 



were also calculated to characterize the spatial dimension 
the polymer chain. Finally we introduced three param- 
eters to measure the degree of helicity within the chain. 
This is accessed in many references by examining the 
Ramachandran angles 4> and ip which are parameters not 
present in our current model. Equivalent to the same 
treatment, the first helicity parameter, Hi, is based on an 
examination of the distance between the i th and (i + 3) th 
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monomers, which is calculated and evaluated as being 
in a helical or non-helical state, as the three connected 
bonds in the segment will be arranged in a very spe- 
cific way when in the helix state. A window of allowed 
distances can then define a three-monomer segment as 
helical, and the average number of helical segments, Hi, 
can be measured. A second method of characterizing the 
degree of local helicity within the polymer involves sum- 
ming the dot products of the adjoining cross products of 
connected bond vectors: 



H 2 = 



1 
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A third measurement is the correlation function G(i) = 
Ui ■ M m id, where u m id corresponds to the vector in Eq. (2) 
for the central monomer, characterizing the correlation 
of the helical ordering along the chain. Discussed below 
is the integrated correlation function, 



H 3 = 
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N-2 
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which represents the averaged correlation in the chain. 
These three parameters yield a global representation of 
the amount of helicity in the polymer. 

Using the reweighting function from preliminary de- 
terministic multicanonical runs, we performed a typi- 
cal production run of 5 x 10 s pivot rotations with mea- 
surement made every 10 steps for rescaled temperatures 
ksT/t ranging from 0.1 to 5.0. Polymer chain lengths 
of 13,19,26, and 39 monomers were considered and the 
results are shown in Fig. 2. 

Displayed in Fig. 2a, the scaled heat capacity curves 
shows a number of interesting features. The first and 
most obvious is the large peak near kT/e = 1.2 corre- 
sponding to the transition to a helix state. The words 
describing what is meant by a helix state must be chosen 
carefully for the reason that there are two distinct helix 
states observed in this study. Figures lb and lc are typ- 
ical configurations just before and after the transition. 
Corresponding to a significant increase in the three he- 
lical parameters measured for N — 19 (Fig. 2c), a C v 
anomaly occurs at this temperature, confirming that the 
polymer is making a transition from a disordered (Fig. 
lb) state to the more ordered helical state (Fig. lc). 
The same sharp increase in C v was also predicted in the 
Zimm-Bragg theory ^[TJ and its variations, and in recent 
MC simulations by Okamoto and Hansmann 0. Zimm 
and Bragg argued that this peak does not represent a 
second order transition in consistent with the known con- 
clusion regarding the impossibility of phase transitions in 
a one dimensional system. However, one can also argue 
that a helix polymer is not strictly a one dimensional 
system. Our peaks in heat capacity curves continuously 
increase as N grow, which show a similar behavior as heat 



capacity curves in a protypical second order phase tran- 
sition, and qualitatively agree with the observation in a 
more realistic polypeptide model for shorter helical loops 
. It would be interesting to further examine a possible 
finite size scaling in the C v peak values determined from 
our simulations. However, currently, the main obstacle to 
achieving this goal is the increasing computational time 
for larger N. 

One of the most interesting features of Fig. 2a is 
the previously unknown second transition occurring in 
the neighborhood of scaled temperature ksT/e = 0.3. 
In a recent study of the phase transitions in isotropic 
homopolymers with specific ground states in a lattice 
model, Zhou et al. |^| have observed a similar transi- 
tion, and have suggested that this is a solid-solid transi- 
tion accompanying the crystallization of their polymers 
into the ground state. In the present case, our wormlikc 
chains have the specifically constructed ground state of 
a perfect helix, and we believe that the second transi- 
tion corresponds to the crystallization of the helix I state 
into a perfect helix II state. Above this transition we 
have observed that the chain is in a highly ordered he- 
lical structure only segmentally as demonstrated by Fig. 
lc. In particular, the end portion of the polymer may 
not necessarily ordered. Even for the ordered segments, 
elongation motion and bending of the ordered segment 
as a whole would make the directional ordering weaker 
compared with the more tightly wound helix II. Thus, as 
the temperature continues to drop further, the few un- 
bonded end segments arrange themselves into a perfect 
helical state (helix II) , and in the mean time, the already 
existing helix segments crystallize to a near ground-state 
structure, shown in Fig. Id. Such a crystallization tran- 
sition has not been seriously modeled by previous theo- 
retical studies, which usually treat the helix I to helix II 
transition as a smooth crossover with no anomalies in the 
heat capacity, as a result of influence by Zimm-Bragg's 
original model fl| -f|jl3l . The three lower curves in Fig. 
2c arc derivatives of the helical parameters, all showing a 
significant change in the slope corresponding to the low- 
temperature peak at the C v curve. To further confirm 
this crystalline transition we have examined individual 
local the helical parameters Ui ■ Ui+i, which also show 
substantial changes before and after the helix I to helix 
II transition. As N becomes large, the peak position in 
C v seems to be to move to higher temperatures, with 
the C v curves becoming flatter. Presently, we are unable 
to confirm these trends due to the large error bars as- 
sociated with the low temperature behavior. Returning 
to Fig. 2a, we would like to point out a third feature 
of interest, which is the emergence of a shoulder in the 
larger polymers near ksT/e = 2, not easily detectable 
from Fig. 2a. This shoulder corresponds to a 9 tem- 
perature collapsing of the polymer chain from a gyrated 
coil ("gaslike") to a randomly collapsed globular state 
("liquidlike"). This collapsing is a direct cause of the 
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FIG. 2. Variation of physical quantities as a function 
of scaled temperature T = fcsT/e. Triangles, diamonds, 
squares, and circles in Figs 2(a) and 2(b) corresponds to 
39mers, 26mers, 19mers, and 13mers respectively. Only ev- 
ery tenth data point was plotted, and more data points are 
represented by the more smooth curves. The error bars of 
the C v curve are almost twice the size of the plotted sym- 
bols, and triple the size of the plotted symbols in the inset. 
Triangles, squares and circles in Fig 2(c) represents the helix 
parameters Hi , H2 and H3 for the case of 19mers, respec- 
tively. The solid, short-dashed, and the dot-dashed curves 
represent dHi/dT, dH^jdT, and dH^/dT respectively. An 
overlay of the C v curve (to the right scale) is also shown by 
the long-dashed curve for reference. 



directionally averaged attractive potential between the 
monomers, creating globular states which have little he- 
lical ordering. In Fig. 2b, the radius of gyration curves 
shows a significant decrease near a transition tempera- 
ture ksT/e = 2, as the signature of a collapsing transi- 
tion in polymers. At lower temperatures, as helix states 
start to form, there is an increase in both R g and i? C nd 
(not shown). Such a turnover is actually a characteri- 
zation of the helix formation, because in a large poly- 
mer, the ordered helix structure will be quite extended 
along the helix axis, forcing the monomers to have an 
average radius of gyration larger than that of a compact 
collapsed state. For shorter helices, the increase of poly- 
mer dimensions at low temperature does not exist, and 
the coil-globular and globular-helix I transitions merge 
to a single transition. Note that neither dimensional pa- 
rameters show appreciable changes at the helix I-helix II 
transition, which is again consistent with our belief that 
the helix I to helix II transition merely involves the or- 
dering of the end segments and tightening of the helix 
segments, with no significant structural rearrangement. 

To summarize, we have described a simple model which 
can be used for examining the coil- helix transitions. To 
furnish a stable perfect helix ground state, one employs 
a strongly directionalized attractive potential. The coil- 
helix transition actually consists of three steps as the 
temperature is lowered. At high temperatures the poly- 
mer has a gaslike coil conformation which collapses into 
a liquidlikc globular conformation with the relative bond 
directions still remain disordered. This state then makes 
a transition to a relatively well ordered "solid" -like helix 
I state, with a dramatic increase in bond direction corre- 
lation. Upon further cooling this conformation crystal- 
lizes to the very well ordered helix II state. It would be 
interesting to confirm that this transition indeed accom- 
panies a C v anomaly experimentally in short polypeptide 
chains which are known to exhibit perfect helical ground 
states. We hope the current model will shed some light on 
analyzing experimental data in complimentary to a tra- 
ditional Zimm-Bragg type physical picture. We believe 
that the current treatment can be modified with little ef- 
fort to produce a variety of ground states, important to 
modeling other interesting transitions in biological sys- 
tems. 
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